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ABSTRACT 

We present new radial velocity measurements from the Bulge Radial Velocity Assay (BRAVA), a 
large scale spectroscopic survey of M-type giants in the Galactic bulge/bar region. The sample of 
~4500 new radial velocities, mostly in the region —10° < / < +10° and b w —6° more than doubles 
the existent published data set. Our new data extend our rotation curve and velocity dispersion profile 
to +20°, which is ^^2.8 kpc from the Galactic Center. The new data confirm the cylindrical rotation 
observed at —6° and —8°, and are an excellent fit to the Shen et al. (2010) N-body bar model. We 
measure the strength of the TiOe molecular band as a first step towards a metallicity ranking of the 
stellar sample, from which we confirm the presence of a vertical abundance gradient. Our survey finds 
no strong evidence of previously unknown kinematic streams. We also publish our complete catalog 
of radial velocities, photometry, TiO band strengths, and spectra, which is available at the IRSA 
archive: http : / / irsa . ipac . caltech . edu/ as well as at UCLA; http : / /brava . astro . ucla . edu/. 

Subject headings: surveys — stars: abundances, distances. Population II — Galaxy: center 



1. INTRODUCTION 

Only for stars in the Milky Way is it currently possible 
to examine both the three dimensional kinematics and 
composition of a central bulge/bar population, offering a 
unique laboratory for the study of galaxy formation and 
evolution. Up until now, there have been very few opti- 
cal spectra and radial velocities of bulge stars, published 
in a catalog form, with the spectra and measurements 
made publicly available. Here we present a catalog of 
the Bulge Radial Velocity Assay (BRAVA) that has a 
sample of ~ 10, 000 M giant stars selected from the red 
giant branch (RGB) of the Two micron All-sky Survey 
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(2MASS; 'Skrutskie etallllOOl). This catalog gives the 
2MASS magnitudes and positions, our measured radial 
velocities, and our spectra. This publication presents our 
complete sample of low resolution spectra for the BRAVA 
survey, which covers approximately —10° < I < +10° 
and —4° < b < —8°. We also present a rich new data set 
at 6 = —6° that complements the first results presented 
in Howard et al. (2008, hereafter Paper I) and the prop- 
erties of which conform to our other latitudinal studies 
at —4° and —8°. We expect that the database will ulti- 
mately be useful in constraining the dynamical model of 
the Galactic bulge, and in placing limits on the detection 
of cold streams from fossil infall events. The 2MASS as- 
trometry is of sufficient quality that proper motions for 
our sample will be measurable in the next few years. 

Although evidence for a barred p otential has long been 
known (e.g. iLiszt fc Burton|[l980l ). the modeling of the 
2 micron s urface br ightness distribution in the bulge as 
a bar (Bli tz fc Sperg cl 1991) established that the stellar 
distribution of the bulge is best modeled as a bar with 
the position angle of its major axis as 20 — 45°. The 
self-consiste nt dynamical model of a rapidly rotating bar 
(|Zhaolll996D was a significant advance, and a number of 
models have followed, however, additional progress re- 
quires constraining new models with a much larger kine- 
matic data set. The question of exactly how the bar 
structure formed - presumably via some kind of secular 
evolution (jKormendv fc KennicuttI [2004D . but including 
rapid star formation, winds, and dissipation - can be 
better addressed with these new large samples. Addi- 
tionally, the relationship between the bulge, inner disk, 
thick disk and halo, may be illuminated further. As the 
coverage and sample sizes analyzed in the bulge increase, 
it becomes possible to place meaningful limits on the in- 
fluence of infalling satellite and other types of merger 
events. The ultimate goal is to move beyond the char- 
acterization of the bar's morphology and orientation to 
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the characterization of the bulge/bar as a stehar system 
with a unique and complex hist ory. 

Early surveys of the bulge (|Nassau fc Blancol [TOSl 
immediately revealed a distinguishing characteristic, 
namely large numbers of M giants compared with the 
giant population in globular clusters and the halo. The 
breakthrough in achieving a physical understanding of 
the M giant population arose from the surveys under- 
taken by Victor and Betty Blanco using the newly com- 
missioned grating/prism on the 4m telescope at CTIO 
(Blanco, McCarthy, & Blanco 1984) - an effort that 
yielded thousands of photographic low resolution classifi- 
cation spectra. These studies provided the raw mat erial 
that fed the first rad ial velocity studv (lMouldlll986D and 
infrared surveys (e.g. lFrogel fc Whitfordll 19871) . The cool 
bulge M giants could be easily identified and exploited 
because of their brightness in the near IR (where extinc- 
tion is low) and by their strong molecular (chiefly tita- 
nium oxide [TiO]) bands. Sharpies, Walker & Cropper 
(1990) used an early fiber spectrograph on the Anglo- 
Australian Telescope to obtain ~ 300 spectra of bulge 
M giants in Baade's Window. However, no study specifi- 
cally exploited the M giants as kinematic probes over the 
entire bulge, even though they could be easily identified 
from their extremely red V — I colors. 

Followi ng the self consistent, rapidly rotating bar 
model of () Zhao! [19961) . two efforts modeled n ewly avail- 
able, substantial data sets. iSevenster et al.l (1999) em- 
ployed the OH/I R star popula tion to constrain a bar 
model, while the .Beaulieu et al.l (|200Q,) study used plan- 
etary nebulae (PNe) as kinematic probes of the bulge 
and fit a range of bar models to the data. Both probes 
offered an advantage in that they are ubiquitous and 
easy to identify over the whole of the bulge. In the 
case of the PNe, the numbers are relatively few (373). 
Unfortunately, PNe are not specific to a given stellar 
population so disk contamination is a possibility, and 
the t otal number in the P Ne sample is modest. How- 
ever, iBeaulieu et al.l ()2000D were able to test both self- 
consistent an d N-body bar models a gainst the PNe kine- 
matic data. ISevenster et al.l ()199 9^ were able to draw 
significant constraints from the OH/IR star kinematics, 
finding a bar-shaped bulge at a roughly 45° bar angle. 
The OH/IR stars are also rare, with numbers similar to 
those of the PNe. 

With the release of the 2MASS database, it became 
possible to easily select M giants over the entirety of the 
Galactic bulge. The 2MASS survey offered an essentially 
unlimited supply of kinematic probes. The Bulge Radial 
Velocity Assay or BRAVA (see also Rich et al. 2007b 
and Paper I) was conceived to exploit the high qual- 
ity, uniform, photometry and astrometry for this data 
set. We select red giants from the K, J— K red giant 
branch that are approximately at the distance to the 
Galactic center, Ro=8 kpc, and are highly likely to be 
bulge/bar members. These stars can thus be selected 
from highly obscured regions of the bulge. The first re- 
sults of the BRAVA project are given additionally in Rich 
et al. (2007b), Paper I, Howard et al. (2009), and Shen 
et al. (2010). The principle result, based on ^4500 stars, 
was the confirmation of cylindrical rotation, and that the 
simple bar/boxy bulge model matches the BRAVA kine- 
matics strikingl y well with no ne ed for a merger-made 
classical bulge (|Shen et al.l I2010D . The sample of red 



giant kinematic probes is now 8585. We also have in 
progress a new version of the self-consistent model (Wang 
et ah 2011). 

As described in lRich et al.l (|2007b( ). we used the Hydra 
spectrograph on the Blanco 4m telescope, set at a central 
wavelength of ■~7900j^[3- Many details of the observation 
and analysis are given in Paper I and will not be repeated 
here. After some experimentation, we ultimately deter- 
mined that even with the prominent, partially overlap- 
ping, molecular band absorption, the near-infrared Ca 
triplet (CaT) offered the best opportunity to obtain ex- 
cellent radial velocity cross correlation. 

While originally conceived as a purely kinematic study, 
the data quality of BRAVA allows us, in principle, to in- 
vestigate the metallicity distribution of the bulge. Only 
with the addition of chemical abundance information will 
we ultimately be able to pin down the formation mecha- 
nisms that laid the kinematic traces found in our survey. 
In particular, the presence of a radial or vertical metal- 
licity gradients in the bulge m ay be indicative formation 
mechanisms jM clcndcz ct all 120081 iZoccali et all I2008I : 
I B abusiaux et al .11 2 1 Q : iJohnson et al. |2011[ ) and warrants 
a well calibrated, large-number sample, such as BRAVA. 
Here we report the coordinates, photometry, TiO molec- 
ular band strengths, and radial velocities of individual 
red giants in the BRAVA fields, as well as the mean ve- 
locities and velocity dispersions for each field. The data 
release also includes all of our actual reduced spectra. 
It is our intention to include proper motion data to the 
compilation, as well. The present paper describes the 
final data set of BRAVA and presents the BRAVA data 
release and website, which includes all available data and 
spectra from our obs erving campaign. We also refine the 
analysis presented in lShen et al.l (|2010) and confirm that 
the Galactic bulge appears to consist of a massive bar un- 
dergoing pure cylindrical rotation and yield strong limits 
on the presence of a classical bulge. 

The structure of the paper is as follows. We present 
the observations and spectroscopy, especially insofar as 
these were changed from the setups used by Paper I and 
iHoward et al.l ()2009l ) in §1. Velocity calibrations and the 
main results are discussed in the next section, §2, where 
we also introduce our public data release. Our molecular 
line index measurements are then introd uced in §3. §4 
shows the new data fit to the iShen et al.l (2010) N-body 
bar model, before concluding the data release in §5. 

2. NEW OBSERVATIONS AND SPECTROSCOPY 

The observations presented here were taken in 2008 
with the Hydra multifiber bench spectrograph at the 
Cassegrain focus of the Blanco 4m telescope at the Cerro 
Tololo Inter-American Observatory (CTI0)E3- Nineteen 
individual bulge fields were observed in July 2008 and 
twenty-three in August 2008. Additionally, three ra- 
dial velocity standard stars (HD 177017, HD 218541 and 
HD 146051), already used in Paper I, were acquired dur- 
ing these runs, as well as observations of the BRAVA 
calibration field at (/, 6) — (6, —4). 

The wavelength range for each BRAVA field is given in our 
data release, as it varies slightly depending on what year the ob- 
servations were taken. 

CTIO is operated by AURA, Inc., under contract to the Na- 
tional Science Foundation. 
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Howard et al. (2008, hereafter Paper I) describe the 
observational setup and sample selection to obtain ra- 
dial velocities of the bulge M giants. Briefly, in 2005, our 
central wavelength was -7600A, but every year it was 
adjusted redward in order to observe the calcium triplet 
as well as the TiO band at ~7050. The specific wave- 
length range for each BRAVA field is listed in our data 
release. For our 2008 observations, the KPGLD grating 
was employed, blazed at 8500A and with 1x2 binning. 
The central wavelength was '^7900A and the effective 
dispersion was 0.88 A pixel"^ with an effective resolving 
power of i? ^ 4200. The spectral range of the new data 
set included all three CaT lines, and for each field, three 
exposures at 600 s each were obtained. These fields typ- 
ically allowed for ^106 fibers to be place on M giants 
and an additional '^20 fibers to be used for sky subtrac- 
tion. The signal-to-noise (S/N) at ^7500 A ranges from 
10-80, with median values of ^ 35-45. The variations in 
S /N arise mainly from transparency at the time the data 
were taken, and the field position in the bulge (extinc- 
tion, crowding). 

A complete listing of the observed fields is shown in 
Table [U and all the observed fields from Paper I as well 
as those presented here are shown in Figure [1] 
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Fig. 1. — All the observed BRAVA fields, overplotted on the 
COBE 2 fj,m image (Launhardt, Zylka, & Mezger 2002). Each 
field is 40' and are color coded to designate the year observed. The 
black circles are fields observed over multiple years. 

The M giants were selected from the 2MASS catalog; 
hence uniform and consistent astrometry and photome- 
try is known for all targets. Figure [5] shows the color- 
magnitude diagram of the targets. Here th e reddening 
values are taken from 'Schlegel et al .l (Il998f). Our tests 
showed that application of the Schle gel et al.l (I1998D val- 
ues led to a consistent decrease in the width of the RGB 
in lower latitude fields, and gave a good overal consis- 
tency between the various fields, which span a wide range 
in Galactic latitude and longitude. Reddening and the 
abundance of target stars decreases as a function of galac- 
tic latitude, so the fields at b=—8° reach to K'^ 10.1 mag 
whereas the fields at b——A° and 6=— 6° go to K~9.5 mag. 
This s election was shown in Paper I and bv lHoward et al.l 
(|2009f ) to largely avoid the red clump in the bulge as well 
as the M giants belonging to the disk. 

The data were reduced and the radial velocities were 
obtained in a similar manner as described in Paper 




(J-K)„ (mag) 



Fig. 2. — The dereddened 2MASS color-magnitude diagram of 
the new 4352 M giants for which radial velocities have been deter- 
mined. Their position on the CMD indicates they are bona fide 
bulge M giants. 

I. After trimming, overscan- and bias correcting the 
spectra, cosmic rays were removed usin g the Laplacian 
edge- detection routine, L. A. Cosmic (van Dokkum et all 
I2001| ). The IRAF task dohydra was used for aperture 
identification and tracing, flat-fielding, wavelength cali- 
bration and sky subtraction. The cross-correlation rou- 
tine, xcsao, in IRAF was then used to cross-correlate the 
spectra against our three radial velocity standard stars, 
obtain the radial velocities and correct them to the he- 
liocentric rest frame. The radial velocity standard stars 
were also used to apply the appropriate zero-point shift 
of the velocities. 

To assess the consistency of our velocity results, a star- 
by-star comparison of stellar velocities was conducted of 
the field at (/, b) = (6, —4) as well as overlap stars from 
the field at {l,b) ~ (0,-6). Figure |3] shows 101 stellar 
velocities in field (/, b) = (6, —4) and the 16 velocities in 
field at {I, b) = (0, —6) that are in common between the 
2006 and 2008 data. As these velocities agree to within 5 
km s^^, (the la dispersion of the difference), this is the 
value adopted as the global, individual stellar velocity 
error. This error is identical to that found for the radial 
velocities in Paper I. This error is also similar to the 
velocity precision obtained with our spectra as reported 
by xcsao. 

The heliocentric velocities are corrected fo r the reflex 
motion of the Sun following iBeaulieu et al.l (f2000i) and 
the velocity dispersion is given by 

N 

^Intrinsic = ^^Ibserved ~ ^[error'^{v^)]/[2{N - 1)], (1) 
i=0 

where cflf^gerved observed velocity dispersion of a 

given field. The second term represents the uncertain- 
ties in the individual stellar radial velocity (vz) mea- 



4 



Kunder et al. 



20 



m 



- 



o 
o 
c\j 
I 

cc 
o 
o 

> 
<] 



-20 



o 2006/2008 
• 2006/2008 





= (6°, 


-4°) 




= (0°, 


-6°) 



o 



0@)o 

o 




o 



o 



o 



o, 



o 



o 



o 



o 



Offset = -2.5 km s-i 
a- 5.1 km 

I I I I I I I I I I I I I L 



9.5 



10 



10.5 



11 



J [mag] 



Fig. 3. — The offset in stellar radial velocities between Paper 
I observations and the 2008 observations for field (Z, b) = (6, —4) 
(open circles) and for field {I, b) = (0, —6) (closed circles). The size 
of the errors in the 2008 velocities as reported by xcsao is shown 
in the top left corner. 

surements, where N is the number of stars in that field. 
Given the large amount of stars in each field, the second 



term is negligible, and hence af^ 



observed ' 



The 



average heliocentric and galactocentric velocities of each 
field is listed in Table [TJ A 4 cr-clipping algorithm is used 
to obtain these values; only one star is affected by this 
clipping. 

The complete BRAVA catalog can be ac- 
cessed at the Infrared Science Archive (IRSA) 
|http : / / irsa . ipac . caltec h . edu/| as well as at 
the University of California, L os Angeles (UCLA) : 
http : //brava . astro . ucla . edu/[ Each BRAVA field 
is designated by its Galactic / and 6, followed by the 
average helocentric and galactocentric velocity, error 
in velocity, dispersion, error in dispersion, as described 
above. Lastly the wavelength range of the spectrum 
is given, as well as the number of stars in each field. 
The catalog includes all data from 2005, 2006, 2007 and 
2008. Following the link beside each field allows the 
measurements of the individual stars in the field to be 
obtained, as well as the reduced spectra in fits format. 
An example of the format for the table of the individual 
radial velocities is shown in Table [H 

2.1. Color /Magnitude Bias 

It has been shown in Paper I and lHoward et ahl (j2009f ) 
that bulge radial velocities summed over minor- and 
major-axis fields consist of a Gaussian distribution with 
no apparent deviation from a normal distribution. The 
data analyzed in Paper I is located physically closer to 
the Galactic center than the new data presented here. As 
the latter is located closer to the bulge/halo boundary, 
a comparison between the data presented here and that 
shown to be from a bona fide RGB bulge sample in Paper 
I can be used to investigate the possibility of color and/or 
magnitude bias in the present sample. The bulge RGB 



stars from the b — —8° and b — —6° fields are summed 
and shown in Figure |4l They yield an apparent Gaus- 
sian distribution with < Vgc > = ~8 ± 3 km s^^ and 
cr = 96 ± 2 km s^i for the b = -8° fields and < Vgc > 
= -10 ± 3 km s^i and cr = 107 ± 3 km s'^ for the 
6 = — 6° fields. Both of these curves are consistent not 
only with each other, but also with those from Paper I, 
which suggests that these bulge samples in fact consist of 
a homogeneous, normally distributed stellar population. 
Any subpopulation is not significant enough to cause a 
departure from the Gauassian distribution. 
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Fig. 4. — Histogram of all bulge RGB star galactocentric veloc- 
ities for the fields at fe=— 8° (bottom) and b = —6° (top). Both 
distributions are normally distributed, with negligible skew and 
kurtosis, which is consistent with a single kinematic population. 



3. TlOe AS METALLICITY INDICATOR 

While our spectral range, in principle, offers a wealth 
of metallicity indicators, the standard calibrations of the 
line strengths of the three prominent Ca II lines at 8498, 
8542^ 8662A onto an [Fe/H] scale (e.g., iRutledge et"all 
119971 ) fail for the cool M-giants that make up the ma- 
jority of the BRAVA sample. Here, the largest prob- 
lem for measuring the CaT is its contamination with the 
strong TiOe band at 8430A, which becomes progressively 
stronger with decreasing Teff. In fact, this band has a 
strong dependence on s tellar effective temperature (e.g. 
IMilone fc BaAuvlll99l . Since we can estimate our Teff- 
scale using the 2MASS infrared colors (Figure IH]) , mea- 
surements of the band strength the n allows to explore th e 
variations with stellar metallicities (jSharples et al.ll990D . 
In particular, at an expected mean of [Fe/H]~ —0.5 dex, 
the TiOe band will be very prominent in the metal rich 
bulge stars. 

As a first step towa rds a comprehensive m etallicity dis- 
tribution, we follow ([Sharpies et all 119901) in extracting 
a TiO-strength index, TiOe, from our spectra, which is 
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defined as a magnitude via 



TiOe = -2.5 log F1/F2 



(2) 



where the fluxes are obtained from straight integration in 
the bands from 8370-8420 A (Fi) and 8440-8490 A {F2). 
The la error on this number, estimated from the vari- 
ance in those spectral bands, is typically 0.05 mag. Val- 
ues of each star's TiO index are provided in the final 
BRAVA catalogue^- Note that this molecular band is 
not measurable in the spectra taken in 2005 due to the 
instrumental set-up. 

Figure [5] compares our TiO measurements for the two 
sets of overlapping spectra, as is already done for the 
velocities in Fig. 3. The la scatter between both sets 
of TiOe is 0.04 mag with a mean difference below 0.01 
mag. This is fully consistent with the uncertainty on the 
index, estimated above from the spectral variance. 
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Fig. 5. — The offset in TiOe between stars observed in botti 2006 
and 2008. Ttie filled circles correspond to stars in field {l,b) = 
(0,-6), and the open circles correspond to stars in field (l,b) = 
(6,-4). 

In Figure [5] the TiOe index is shown as a function of 
infrared color. The distribution of the molecular index 
in our spectra shows a narrow band around zero, i.e., 
there is no discernible TiO depression, but a well-defined 
pseudo-continuum in the warmer stars. This holds for 
a broader range in J— K between 0.7 and 1.1 mag. On 
the other hand, for stars redder than J— K>1.0 (corre- 
sponding to Toff<3700 K) there is a notable presence of 
progressively strong bands, leading to a linear rise in the 
TiO index. 



3.1. Phc 



models 



A dispersion in TiOe at fixed (J — K) can be due to 
both differential reddening and metallicity spread. To 
illustrate this, we consulted a set of synthetic spectra 
based on the Phoenix stellar atmosphere models. These 

In a future work (Koch et al. in prep.) wc will explore the 
reliable calibration of the CaT line strengths onto metallicity for 
a subsample of the stars with warmer temperatures, in which the 
CaT is yet unaffected by TiO absorption. Our measured equivalent 
widths will then be added to this BRAVA release. 



Fig. 6. — TiOe index (after Sharpies et al. 1990) versus infrared 
color. The solid lines in the bottom panel are based on a-enhanced 
Phoenix models for metallicities, [M/H] between —2.5 and +0.5 dex 
(right to left) in steps of 0.5, adopting log(/=l. A typical error bar 
is indicated in the top left corner. 

offer a large grid of stellar atmospheric parameters and 
a spectral library in high resolutiorP^ (Hauschildt, Al- 
lard & Baron 1999, Hauschildt et al. 2003). In particu- 
lar, its meticulous treatment of cool stellar atmospheres 
and its vast input sets of molecular opacities makes this 
synthetic library ideally suited for comparison with the 
BRAVA data. Our chosen model grid consisted of T^g 
ranging from 3200-4500 K in steps of 100 K, logg = 0, 
0.5, 1.0, 1.5, metallicities [M/H] = -2.5 to +0.5 in incre- 
ments of 0.5 dex, and a-enhancement of -1-0.4 dex as well 
as solar-scaled abundances, i.e., [a/Fe]=0 dex. This gen- 
erously comprises the range of parameters spanned by 
the BRAVA late-type giant sample and is also consistent 
with the range in spectral types targete d in abundance 
studies of comparab le stellar s amples (IRich fc Origlial 
I200I IFulbright et"alll2006.: Rich etani2007aD . The swT 
thetic spectra were then smoothed with a Gaussian ker- 
nel to mimic the resolution of the BfiAVA spectra. Fi- 
nally, all the relevant quantities such as the TiOe index 
were derived from these convolved spectra in the identi- 
cal manner as for the observed data. 

As Figure [6] indicates, the bulk of the target stars 
falls near [M/H] of —0.5 dex, which is consistent with 
the metallicity distribution functions in the literature 



(e.g. , RanEwala & Williams"200l; iDe Propris et al.ll201]] 
[Johnson et al. 2011). Moreover, stars are found to cover 
a broad range from super-solar down to low metallicities 
around —1.5 dex, with a hint of a very metal poor ex- 
tension towards —2 dex. A detailed, quantitative study 
of the abundance properties of the BfiAVA sample will 
be presented in a forthcoming paper. 

A number of points on the blue side of Figure |6] seem 
inconsistent with th e Phoenix models. One reason for 
this may be that the lSchlegel et al.l (|1998l) reddening val- 
ues overestimate the reddening, and these points should 

ftp:/ /phoenix. hs .uni-heunburg . de/GAIA.HighRes/SpectraT' 

Note, however, that such observations strictly derive MDFs 
based on iron abundances, [Fe/If], and a quantitative comparison 
with the models' global metallicity may not be straighforward. 
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then come from the most reddened fields. We should ex- 
plore whic h fields accoun t for t hese points. When com- 
paring the lSchleeel eFall (flMl E(B-V) values f or 330 
BRAVA stars located within the iGonzalez et all (|2011l ) 
reddening m ap, a ^0.15 mag o ffset in E(B — V) is seen, 
with the iSchlegel et al.l ()1998f ) values being systemati- 
cally larg er. This translates to ~0.08 mag in (J — K); 
hence the IGonzalez et al.l ()201lD reddening values would 
move the (J — K)o values in Figure|6]to the red by roughly 
the spacing between the Phoenix curves, and the number 
of points on the blue side of the -1-0.5 dex Phoenix model 
would be decreased. Another possible contribution of un- 
certainties in the TiO strength is the uncertainty in the 
2MASS photometry, which is roughly 0.03 to 0.1 mag in 
( J — K) . The re has also been the long standing problem 
in the bulge (jFrogel fc W hitford 1987) that the bulge gi- 
ants appear to have TiO that is much too high for their 
(J— K). Further surveys such as the APO Galactic Evolu- 
tion Experiment (APOGEE) will obtain high-resolution, 
high S/N infrared spectroscopy of red giants stars across 
the Galactic bulge, and may shed light on this issue. 

3.2. Metallicity gradients 

Recent high resolution spectroscopic studies of the 
bulge find a vertical abundance gradient in the sense 
of lower metallicities towards h igher latitudes (e.g^ 
Zoccah et all 120081: iJohnson et "all 120111: [Gonzalez et atl 
201 It) . Based on the large BRAVA sample, this issue 
can also be investigated (although, here, we only do it 
qualitatively). Figure [7] shows the TiO metallicity indi- 
cator separated by Galactic latitude. The Galactic lon- 
gitude of these fields is /=0°, and the Phoenix model for 
a metallicity of —0.5 dex is over-plotted to guide the eye. 
It becomes immediately obvious that there is in fact a 
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Fig. 7. — Same as Fig. (6] but separated by the Galactic latitudes 
covered by BRAVA. The gray-shaded region indicates the distribu- 
tion in the 6 = — 4° field, and the Phoenix model for a metallicity 
of —0.5 dex is over-plotted. 

systematic difference between the fields with a tendency 
for the lower latitude fields to have higher metallicities 
- a trend consistent with recent suggestions of metallic- 
ity gradients across the bulge derived from the K giants 
(e.g., Zoccah et al. 2008; Johnson et al. 2011). This 



is also consistent with results of the M giant population 
from Frogel, Tiede & Kuchinski 1999; they use IR CMDs 
to show there is a gradient in mean metallicity along the 
minor axis in the range 0° > 6 > —10°. Moreover, we find 
the most prominent trend concerns the increase in metal- 
licity spread as one approaches the inner bulge regions. 
However, a detailed analysis on the uncertainties in the 
reddening and photometry in the inner bulge regions as 
compared to those in the outer regions is necessary, espe- 
cially as th is is not seen in the high resolution data (i.e. 
IZoccah et al. 2008; Johnson et al. 2011). 

Finally, we show in Fig. [5] the entire sample segregated 
by stellar (K-band) magnitude. Amongst all of the mag- 
nitude bins the color-TiO distributions are broadly con- 
sistent with each other, bolstering our findings from Sec- 
tion 2.1 that argue against any color/magnitude bias in 
the BRAVA sample - also with regards to metallicity. 
The only apparent exception is the faintest bin, which, 
however, exclusively contains stars at the highest latitude 
(Figure[2]). Thus, the metallicity gradient permeates also 
in this (magnitude-) subsani ple. 

If the Schleg el et al ] (1 19 98f) extinction were higher than 
that of Gonzale z et al.l (|201lT). it could in principle cause 
our derived ( J — K)o to be too blue in the fields of greater 
ex tinction. As d is cusse d above, the lower E(B — V) 
of IGonzalez et al.l ()201lD have an 0.08 mag effect on 
(J — K)o. However, this is unlikely to mimic a metallicity 
gradient. Small-scale variations in reddening are partic- 
ularly strong in the inner regions, but these regions are 
mostly constrained to |6|<4° (see IGonzalez et 311120111 
their Fig. 6 and Fig. 7). The metallicity gradient in the 
BRAVA data extends all the way to 5=— 1.5°, but is 
also prominent when considering the &=— 8° to 6=— 4° 
fields. Further, we find a similar scatter ('^0.1 mag) 
and a similar o ffset ( ^0.1 5 mag) in AE(B — V) b etween 
ISchlegel et all (| 19981 ) and IGonzalez etaD (|2011f ) of our 
5=— 4° and 6=— 6° stars. Yet these fields still show a 
metallicity gradient. This strongly suggests that redden- 
ing uncertainties can not fully explain away the slight 
metallicity gradient we are seeing. 
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4. OUTLIERS 
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Radial velocity surveys of the Milky Way halo have 
the potential to identify rare velocity outliers that result 
from dynamical processes throughout the Galaxy. The 
most obvious outlier in our sample is 2MASS 17464606- 
3937523, with a heliocentric radial velocity of 447 km/s. 
Its spectrum is shown in Figure [HI The S/N is high, 
and to first order it looks metal poor. Unfortunately, 
at a,6) = (-9.0433,-5.7 386), it is not in the OGLE-II 
proper motion catalog (jSumi et al.ll2004[ ). Hyperveloc- 
ity st ars (HPVs) were first discovered by IBrown et al.l 
(|2005D . and are generally B-type stars moving 2-3 times 
the Galactic escape velocity. At a heliocentric velocity 
of ^450 km/s, the star may not qualify as a true hy- 
pervelocity star, although at r=50 to 100 kpc, unbound 
stars have such velocities (|Kenvon et al.|[2008[ ). A pre- 
hminary [Fe/H] based on the CaT is — 0.86±0.05 dex, 
with an error purely based on the continuum variance 
and magnitude errors, not accounting for uncertainties 
in the calibration coefficients. This is not an unusual 
[Fe/H] value for a bulge star. 
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Fig. 9. — Wavelength calibrated spectrum of a velocity outlier 
in our sample, with a Heliocentric velocity of 447 km/s. The CaT 
lines are labeled in the spectrum for reference. 



5. RESULTS AND ANALYSIS 

Radial velocities, especially in minor axis fields, can 
constrain bulge models. Figure [TU] shows the disper- 
sion profile and rotation curve for the minor-axis strips 
at I = -6°, I = 6°, a nd 1 = 0°. The predictions 
of the iShen et al.l (|2010D models is over-plotted. The 
Shen model uses a cylindrical particle-mesh code to con- 
struct an TV-body model to the BRAVA data. The 
initial parameters adopted are an unbarred disk and a 
rigid pseudo-isothermal halo potential (see Shen et al.l 
120101 for details). In this model, the Milky Way self- 
consistently develops a bar, which buckles and thickens in 
the vertical directions. Hence, there is no classical bulge 
component, and the best-fit model predicts a bar half- 
length of ~4 kpc, extending 20° from the Sun-Galactic 



center line. The additional data presented here allows 
the Shen model t o be compa r ed to two more minor axis 
fields. As seen in lShen et al.l ()2010[ ). the data agree well 
with the Shen model for both the mean velocity and ve- 
locity dispersion observations. 
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Fig. 10. — Velocity dispersion profile (bottom) and rotation curve 
(top) for the /=— 6°, +6° and 0° strips. The open symbols indicate 
data already published and th e filled sym b ols in dicate the data 
presented here. The model of IShen et al.l II2010I ) is overplotted, 
with good agreement to the observations. 

Figure [TT] shows the dispersion profile and rotation 
curve for the major-axis strips at 6 = —4°, b = —6°, and 
b = —8°. The b = —4° data comes from Paper I nd half of 
the b = —8° data comes from Howard et al. (2009). The 
Shen mo del matches t he ob servations well, validating the 
results of IShen et al.l (|2010l ) . It is also worth noting that 
the rotation curves at 6 = —4°, —6°, and —8° are remark- 
ably similar, suggesting that the Galaxy's bulge rotates 
cylindrically. This had been suggested bv iHoward et al.l 
(2009), but here the observational data are more plenti- 
ful and include additional fields at b = —6°. Thus our 
bolstering of the cylindrical rotation signal is further ev- 
idence that the bulge is an edge-on bar, as predicted by 
the Shen model. 

The velocity dispersion drops as one moves further 
in latitude from the Galactic disk plane. A hotter 
model contributes to a large velocity dispersion, and the 
nucleus-spheroid component of the model is hot. How- 
ever, the BRAVA data show remarkable agreement with 
the disk/bar component of the model, with a relatively 
flat dispersion profile at ~70 km contrasting with the 
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Fig. 11. — Velocity dispersion profile (bottom) and rotation curve 
(top) for the 6=-4°, -6° and -8° strips. The filled symbols in- 
dicate data already published and the open symbols indicate the 
data presented here. 

spheroid dispersion, which is predicted to be at ~120 km 
s~^. The turnover seen in the rotation curve at 6 = —4° 
is not evident a.t b = —6° or 6 = —8°. 

Additionally, in Figure [121 the BRAVA fields extend- 
ing past the main body of the bulge at />10° are plotted. 
These observations lie outside the high surface brightness 
boundary of the COBE bulge, and can be used to exam- 
ine the extent of the bar/bulge and to look for signs of 
disk contamination. The full le ngth of the bar is thought 
to be on the order of 3-4 kpc (iHammerslev et al.l 120001 : 
iBissantz fc GerhardI l200l iCabrera-Lavers et al.l l2007f ). 
and since the bar's pattern speed is relatively rapid, it 
effectively "controls" the region that it lives in. At the 
distance of the Galactic Genter (7.9 kpc), 4 kpc at the bar 
angle of ^^20 degrees corresponds to ±1^16°. Indeed, the 
fields at l>15° show a drop in velocity dispersion as well 
as a higher mean velocity, which is due presumably to 
the presence of the inner disk component in these fields. 
Figure [12] also shows the observ ations from a sample of 
373 PNes ([Beaulieu et al.ll2000D : the BRAVA data con- 
firms the rotation seen by the PNes at \l\ > 12°, as well 
as a drop in the velo city dispersion i n the se fields. The 
region covered in the lBeaulieu et aTl (2000) data is —20° 
< / < 20° and -5° > 6 > -10°. Observations at |/| > 
12° place important constraints for Galactic bulge the- 
oretical models and predictions, such as the bar angle 
(Martinez- Valpuesta, private communication). 

The galactocentric velocity distributions of the new 
bulge fields are now presented. As with Paper I, a 



Fig. 12. — Velocity dispersion profile (bottom) and rotation curve 
(top) for all the fe=— 4° fields. The filled symbols indicate data al- 
ready published and the open symbols indicate the data presented 
here. Over-plotted is data from bulge PNes from Beaulicu ct al] 
II2000I) . as well as the predictions of the lShen et al.. (.2010.^ model. 
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bin width of 25 km s^^ is used, driven by the approxi- 
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Fig. 13. — Bulge line-of-sight velocity distributions continued 
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Fig. 13. — Bulge line-of-sight velocity distributions continued 
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Fig. 13. — Bulge line-of-sight velocity distributions continued 

mate dispersion expected for cold components like a dis- 
solving cluster or dwarf gal axy (~10 km/s). Simula- 
tions by iReitzel et al.l ()2007[ ) show that with the num- 
ber of stars we haye in each field (~100), deviations 
from a normal distribution are not unexpected in random 
draws. However, fields with the largest deviations from 
a normal distribution may aid in the selection of follow- 
up observations. To test for normality, the Shapiro- 
Wilk test (Shapiro & Wilk 1965; D'Agostino 1986) and 
the A nderson-Darling test (jPress. Flannerv fc Teukolskvl 
Il986f ) are employed. These normality tests are two of the 
most powerful tests for deviations from normality and are 
especially useful for sample sizes that are relatively small. 
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Fig. 13. — Bulge line-of-sight velocity distributions continued 

i.e., ■~100 particles. For example, the Shap iro- Wilk nor- 
mality test has successfully been used by iKeplev et al.l 
p007i ) to detect known streams in one component of ve- 
locity. For both of these tests, the smaller the p-value, 
the less likely it is that the data come from a normal dis- 
tribution, and a p- value < 0.05 rejects the supposition of 
normality. The results of these tests are given in Table [Sj 
The tested BRAVA field is first listed in Table [3] followed 
by the number of stars in the field, the skewness and kur- 
tosis, and its associated p-value from the Shapiro- Wilk 
test and from the Anderson-Darling test. 

In general, the skewness and kurtosis values are all 
within one or two sigma, the exception being fields {l,b) 
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Fig. 13. — Bulge line-of-sight velocity distributions continued 
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Fig. 13. — Bulge line-of-sight velocity distributions continued 

= (-1,-8), (-4,-8), (-6,-8), (8,-8) and (9,-8). Upon 
closer examination, it is apparent that except for (^,6) = 
(-6,-8) which is discussed in more detail below, these 
fields each contain a star on the tail end of the distri- 
bution which greatly influences the skewness and kur- 
tosis values. These fields also all fail the Shapiro- Wilk 
normality test, which gives more weight to the tails of 
the distribution than the Anderson-Darling test. All of 
these fields are located at 6=— 8°, and it is likely that we 
are seeing some signatures of the bulge / disk boundary in 
these fields. 

There are only two BfiAVA fields that fail both the 
Shapiro- Wilk and the Anderson-Darling test. These 



fields are located at (^,5) = (0,-3.5) and (Z,6) = (-6,-8). 
The former has a skewness and kurtosis consistent with a 
Gaussian distribution, but p-values that reject the sup- 
position of normality. Removing the stars that lie in 
the tails of the distribution does not cause the p-value 
to increase and hence indicate normality. There is a 
strong peak in the velocity distribution at around —30 
km/s and also not many stars at ~100 km/s, which 
are the likely features causing the Shapiro- Wilk and 
Anderson-Darling tests to reject normality. Interestingly, 
the possible signature o f a dis rupted satellite reported by 
iRangwala fc Williams! (pOOOl ) at (Z,6) = (-h5.5,-3.5), has 
a velocity distribution excess around —35 km s^^. How- 
ever, we find no correlation at (i,6) = (0,-3.5) with the 
stars that have velocities around —30 km/s and their TiO 
values. Our sample size is small, and follow-up observa- 
tions of this field would be particularly interesting. 

The field at (/,6) = (-6,-8), on the other hand, has 
the largest skewness and kurtosis of any BRAVA fields 
(1.38 ± 0.25 and 4.91 ± 0.51, respectively). Further, it 
has three stars with heliocentric velocities, Vhc, greater 
than 300 km s~^; the radial velocity dispersion of these 
three stars is 16 ± 7 km/s. From our complete sample of 
8585 stars, only 11 stars have velocities that are above 
300 km/s, and so it is striking that three of these stars 
are located in this field. Again it would be especially 
interesting to obtain follow-up observations of the stars 
in this field to see if more stars with such large velocities 
are found. 

Figure [14] shows all the 8585 radial velocities obtained 
from the BRAVA dataset co- added (bottom panel). The 
mean is 2 ± 1 km s~^ with a cr = 107 ± 1 km s"""^. The 
skewness is negligible (0.03 ± 0.03) and the kurtosis of 
—0.17 ± 0.05 implies a slightly platykurtic (flattened) 
distribution. Both of these results are in agreement with 
the results found by Paper I. Further, the lack of sig- 
nificant skewness and a small value of kurtosis is con- 
sistent with our argument that our Bulge radial velocity 
distribution is not largely contaminated by either cold 
components (disk) or hot components (halo). 

A Shapiro- Wilk normality test fails for this co-added 
dataset, as the sample is too large; with 8585 stars, there 
are many radial velocities with very similar values. An 
Anderson-Darling normality test gives a p-value of 0.59 
and a Kolmogorov-Smirnov test, which works best with 
large datasets, gives a p-value of 0.10. These results sug- 
gest no significant deviations from a Gaussian distribu- 
tion. 

The top panel of Figure [U shows the radial velocities 
shifted onto zero and then co- added. The mean is — 3 ± 1 
km s^^ with a cr = 95 ± 1 km s~^. Again, the distribution 
is Gaussian, with no apparent deviation from a normal 
distribution. It is noteworthy that there is only one star 
in our sample with a velocity greater than ±4 tj, which 
is discussed in §4. 

Figure [15] shows the longitude-velocity plot for the 
three BRAVA major axis strips. For latitudes closest 
to the plane (6=— 4° and 6=— 6°) there is no evidence of 
a cold, disk component in our sample which would man- 
ifest itself as a linear trend. However, further from the 
plane (5=— 8°, ~1.2 kpc from the Galactic plane), the 
"S" shape is not as apparent, suggesting that the disk 
does contribute to the BRAY A sample in this regime. 
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Indeed as compared to the 6=— 4° and b=—6° fields, the 
8 ° fields are also those that in general have skew- 
ness, kurtosis and p- values that are least consistent with 
a normal distribution (see Table [3]) . Recent studies of 
planetary nebulae in t he inner MW suggest a bulge-disk 
interface at ~1.5 kpc (jCavichia. Costa fc Macicl 20111 ). 
consistent with our findings. 

Figure [16] shows the latitude-velocity plot for three 
BRAVA minor axis strips. As expected a linear trend 
is seen; neither cold components nor indications of a hot, 
non-rotating population is seen. However we see "rota- 
tion" in the sense of solid body rotation. 
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-200 
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Fig. 14. — The co-added galactocentric radial velocities of the 
BRAVA stars. In the top panel, each BKAVA field is shifted to 
zero before co-addition. A bin width of 15 km s~^ is used for both 
histograms. Neither distributions show any signs of deviation from 
a Gaussian. 



5.1. A Milky Way Pseudobulge? 

A vertical metallicity gradient i n the bulge has been 
reported in several recent studies (Melendcz et al. 20081 
Zoccali et al. 2008; Babusiau x ct al. 2010; Johnson et aJ] 
201 in . and it has often been used as an argum ent against 
the boxy bulge/bar model of [SheneFaDdloTS)- It is true 
that the collective bar buckling happens to essentially the 
whole disk that is in place at the time of the buckling in- 
stability, but a vertical metallicity gradient could still be 
consistent with the Shen bar/pseudobulge model. One 
possibility is that some of the vertical thickening could be 
prod uced by resonaii t heating of stars that scatter off th e 
bar (|Pfennigeil[l98l [l985l: IPfenniger fc Normfml [l990l) . 
If the most metal-poor stars are also the oldest stars, then 
they have been scattered for the longest time and now 
reach the greatest heights, hence givi ng rise to a verti - 
cal metallicity gradient. Secondly, the lShen et al.l ()2010l ) 
model still allows a relatively small merger-built classical 
bulge (about 10% of the disk mass), the mixing of two 
bulge populations could conceivably produce a vertical 
metallicity gradient (also see the Besancon model, Robin 
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Fig. 15. — Longitude- Velocity (Z„) plot for the entire bulge sam- 
ple at b=~4° (top), fe=-6° (middle) and b=-8° (bottom). The 
plot is smoothed to 1° in longitude and 10 km s~l in galacto- 
centric velocity. This figure shows the cylindrical rotation trend 
very clearly. Notice the lack of any prominent "cold" features that 
would indicate a possible stream detection across multiple fields. 

et al. 2011). A mixing of the bar/pseudobulge with the 
inner thick disk was also proposed to explain t he vertical 
metalicity gradient tB ekki fc Tsuiimotoll201l[ ). 

That the Milky Way has a pseudobulge is consistent 
with the majori ty of the bright galaxies in the local 
11 M pc volume (|Kormendv et al.ll2010l : iFisher fc Drorvl 
120111) . especially for those galaxies with a mass similar 
to that of the Milky Way. Recently an in-depth study of 
NGC4565 has revealed that this galaxy also contains a 
pseudobulge and no hint of a m erger-built classical bulge 
(jKormendv fc BarentintjMol) . Like the MW, NGC 4565 
has a peanut shaped bulge wh ich rotates cylindrically 
(jKormendv fc Illingworthlll982D an d has a strong abun- 
dance gradient ([Proctor et al.ll2000[ ). Hence other galax- 
ies besides the MW with an abundance gradient but no 
(or very little) classical bulge are known to exist. Evolu- 
tion of the vertical metallicity gradient however has not 
been extensiv ely investigated by theore tical studies, but 
the model bv iBekki fc Tsuiimotol ()2011[ ) is a step forward 
in understanding why an abundance gradient is observed 
in the bulge even if the bar evolved secularly. 

6. CONCLUSIONS 

We have presented the final data set of the 
BRAVA survey and provided a website for data 
access, which is available at the IRSA archive: 
http : / / irsa . ipac . caltech . edu/ as well as at UCLA: 
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Fig. 16. — Latitude- Velocity plot for the entire bulge sample 
at /=-6° (top), 1=0° (middle) and 6=+6° (bottom). The sample 
is smoothed to 1° in latitude. 



http://brava.astro.ucla.edu/. The observations at 
the 6° latitude fields as well as more observations at 
the 5=— 8° latitude fields confirms previous suggestions 



that the Galactic bulge rotates cylindrically, a s do boxy 
bulges of other galaxies (e.g., iKormendv fc Il lingwortb] 
1982; 'Jar vislUOOOl: [Sha^ 119931: IKormendv fcKennicutl 
2004; How ard et al.ll2009( ). T he data in this paper double 
the number of stars shown inlHoward et al.l (|2008'. '2OO90 
and used in lShen et al.l (|2010l ). Analysis of these data in 
the same fashion shows that the Galactic bulge consists 
of a single massive bar formed by secular evolution. The 
existence of such a massive bar and no classical bulge 
implies that our Galaxy has not suffered any large merg- 
ers (enough to form a classical bulge) since the epoch 
at which the disk formed. Based on TiOe index mea- 
surements we qualitatively see the signatures of a verti- 
cal metallicity gradient - an attribute that is in contrast 
to our kinematic characterization of the bulge as a pure 
pseudobulge. A detailed analysis of the BRAVA metallic- 
ity scale and its spatial distribution is clearly warranted 
and underway. 
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TABLE 1 

BRAVA Rotation and Dispersion Results 



Date of Obs 


R.A. 


Decl. 


gal I 


gal b 


Target 


Reliable 


< Vhc > 


< Vgc > 


err<v> 


(T (km 




(VT) 


(J2000.n) 


(J2nnn.n) 




f<lrg) 


Fibors in 

(.'onlig 


A'olnril irs 


(km s-M 


(l«n s-i) 


(km s-i) 




(kjlJ 

') 


2008 Jul 08 


17:57:05.75 


-40:43:23.4 


-9.0 


-8.0 


108 


104 


-56.10 


-84.26 


8.32 


84.86 


5.88 


2008 Jul 08 


18:02:09.79 


-38:59:18.5 


-7.0 


-8.0 


102 


102 


-34.14 


-54.39 


7.49 


75.60 


5.29 


2008 Jul 08 


18:14:40.35 


-25:42:48.7 


6.0 


-4.0 


107 


106 


43.38 


76.03 


8.56 


88.12 


6.05 


2008 Jul 08 


18:36:28.77 


-24:54:15.4 


9.0 


-8.0 


109 


108 


11.27 


54.85 


6.95 


72.20 


4.91 


2008 Jul 09 


18:06:51.74 


-37:15:44.5 


-5.0 


-8.0 


110 


108 


-43.43 


-55.69 


7.76 


80.62 


5.49 


2008 Jul 09 


18:11:25.87 


-35:30:15.1 


-3.0 


-8.0 


115 


114 


-36.62 


-40.84 


7.89 


84.27 


5.58 


2008 Jul 09 


18:22:18.76 


-31:05:47.0 


2.0 


-8.0 


105 


103 


-8.65 


7.16 


7.61 


77.21 


5.38 


2008 Jul 09 


18:24:32.95 


-30:12:46.9 


3.0 


-8.0 


109 


108 


-12.17 


7.74 


8.70 


90.39 


6.15 


2008 Jul 09 


18:28:19.59 


-28:29:36.9 


5.0 


-8.0 


106 


102 


35.93 


63.48 


9.04 


91.34 


6.39 


2008 Jul 10 


17:41:55.43 


-36:03:35.9 


-6.5 


-3.1 


65 


63 


-40.86 


-58.03 


11.46 


90.96 


8.10 


2008 Jul 10 


17:50:54.07 


-37:02:24.9 


-6.5 


-5.1 


89 


79 


-33.10 


-50.23 


9.12 


81.09 


6.45 


2008 Jul 10 


18:10:07.40 


-25:39:30.9 


5.5 


-3.1 


88 


86 


26.54 


57.80 


8.95 


83.02 


6.33 


2008 Jul 10 


18:32:36.91 


-26:39:15.7 


7.0 


-8.0 


104 


102 


11.58 


47.41 


7.55 


76.21 


5.34 


2008 Jul 11 


17:59:22.08 


-37:39:25.0 


-6.0 


-6.9 


102 


97 


-45.49 


-61.92 


8.10 


79.79 


5.73 


2008 Jul 11 


18:17:49.98 


-26:10:14.0 


6.0 


-4.8 


107 


106 


13.82 


45.92 


7.83 


80.56 


5.53 


2008 Jul 11 


18:22:15.60 


-26:41:40.5 


6.0 


-6.0 


104 


102 


16.35 


48.22 


8.32 


84.03 


5.88 


2008 Jul 11 


18:26:03.24 


-27:07:50.2 


6.0 


-6.9 


99 


98 


14.94 


46.70 


8.46 


83.78 


5.98 


2008 Jul 11 


18:29:57.01 


-23:06:19.2 


10.0 


-5.9 


93 


88 


20.85 


68.56 


8.72 


81.54 


6.17 


2008 Jul 12 


17:55:51.00 


-37:05:52.7 


-6.0 


-6.0 


102 


96 


-51.06 


-66.83 


8.70 


85.28 


6.15 


2008 Aug 18 


17:45:20.95 


-40:34:09.7 


-10.0 


-6.0 


104 


103 


-53.70 


-85.72 


7.10 


72.09 


5.02 


2008 Aug 19 


17:48:04.97 


-39:41:08.2 


-9.0 


-6.0 


105 


103" 


-47.46 


-75.29 


8.21 


82.91 


5.81 


2008 Aug 19 


17:50:44.04 


-38:51:51.5 


-8.0 


-6.0 


106 


105 


-72.36 


-96.35 


7.27 


74.54 


5.14 


2008 Aug 19 


17:53:05.49 


-38:00:47.4 


-7.0 


-6.0 


111 


103 


-42.71 


-62.86 


9.71 


98.50 


6.86 


2008 Aug 19 


18:26:34.99 


-24:52:48.8 


8.0 


-6.0 


109 


107 


28.44 


68.56 


8.21 


84.88 


5.80 


2008 Aug 20 


17:58:08.76 


-36:17:31.2 


-5.0 


-6.0 


109 


108 


-48.01 


-60.06 


8.43 


87.60 


5.96 


2008 Aug 20 


18:28:33.61 


-23:59:55.1 


9.0 


-6.0 


111 


111 


18.41 


62.45 


7.93 


83.55 


5.61 


2008 Aug 20 


18:30:21.68 


-23:06:39.7 


10.0 


-6.0 


105 


104 


12.49 


60.42 


8.61 


87.81 


6.09 


2008 Aug 21 


18:02:53.91 


-34:31:17.8 


-3.0 


-6.0 


111 


109 


-25.85 


-29.81 


8.33 


86.94 


5.89 


2008 Aug 21 


18:07:17.57 


-32:47:21.0 


-1.0 


-6.0 


111 


110 


-11.74 


-7.74 


9.26 


97.15 


6.55 


2008 Aug 21 


18:11:48.68 


-31:01:50.7 


1.0 


-6.0 


114 


112 


-21.54 


-9.42 


9.16 


96.93 


6.48 


2008 Aug 21 


18:16:10.33 


-29:19:15.1 


3.0 


-6.0 


115 


114 


10.17 


30.18 


8.98 


95.91 


6.35 


2008 Aug 22 


18:00:21.58 


-35:24:53.8 


-4.0 


-6.0 


114 


112 


-43.53 


-51.61 


8.35 


88.41 


5.91 


2008 Aug 22 


18:05:10.43 


-33:40:05.1 


-2.0 


-6.0 


112 


109 


-40.98 


-40.96 


8.55 


89.29 


6.05 


2008 Aug 22 


18:20:20.14 


-27:33:33.3 


5.0 


-6.0 


111 


109 


25.20 


53.17 


8.46 


88.38 


5.99 


2008 Aug 23 


18:09:46.09 


-31:56:01.1 


0.0 


-6.0 


110 


109 


-8.53 


-0.42 


8.70 


90.82 


6.15 


2008 Aug 24 


18:13:54.50 


-30:07:05.2 


2.0 


-6.0 


112 


110 


4.04 


20.30 


8.31 


87.11 


5.87 


2008 Aug 24 


18:18:25.27 


-28:25:38.5 


4.0 


-6.0 


116 


113 


4.99 


29.13 


8.18 


87.00 


5.79 


2008 Aug 24 


18:24:29.51 


-25:47:41.7 


7.0 


-6.0 


111 


111 


19.66 


55.65 


8.36 


88.10 


5.91 


2008 Aug 24 


18:26:51.06 


-20:27:05.7 


12.0 


-4.0 


105 


104 


29.39 


85.59 


7.84 


79.94 


5.54 


2008 Aug 25 


18:30:38.75 


-18:36:13.3 


14.0 


-4.0 


110 


110 


34.63 


98.83 


7.65 


80.26 


5.41 


2008 Aug 25 


18:38:04.33 


-15:04:28.9 


18.0 


-4.0 


106 


103 


28.39 


107.73 


6.09 


61.77 


4.30 


2008 Aug 25 


18:45:38.33 


-11:32:11.6 


22.0 


-4.0 


103 


101 


37.81 


132.19 


6.36 


63.91 


4.50 



"Number of velocities before the cr-clipping algorithm removed one star in the calculation of the mean velocity and 
dispersion. 
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TABLE 2 

Data For Individual Stars in BRAVA Fields 



Aperature 


gal ( 
(deg) 


gal 6 R.A. 
^degj t^JZUOU.UJ 


Decl. 
^J.^OUO.OJ 


J 

(mag) 


H 

(mag) 


K 

(mag) 


V„c 
(km 


E(B-V) 


Jo 

(mag) 


Ho 

(mag) 


Ko 

(mag) 


TiO(mag) Spectrum 
(fits) 


2MASS 
ID 


1 


9.0256 


-7.8183 278.9434 


-24.7927 


10.640 


9.940 


9.724 


-56.4 


0.405 


10.275 


9.707 


9.575 


-0.0345 


Spectrum 


18354640- 
2447338 


3 


8.6962 


-8.0373 279.0005 


-25.1823 


11.116 


10.381 


10.136 


101.8 


0.418 


10.739 


10.140 


9.983 


-0.0390 


Spectrum 


18360011- 
2510562 


4 


9.0018 


-7.8323 278.9456 


-24.8201 


10.468 


9.625 


9.298 


-84.1 


0.403 


10.105 


9.393 


9.150 


0.0055 


Spectrum 


18354694- 
2449123 


7 


9.2062 


-7.7561 278.9694 


-24.6047 


10.643 


9.843 


9.590 


-40.8 


0.382 


10.298 


9.623 


9.450 


-0.0336 


Spectrum 


18355264- 
2436170 


8 


8.9748 


-7.9853 279.0837 


-24.9118 


10.018 


9.043 


8.722 


115.9 


0.418 


9.641 


8.802 


8.569 


0.0895 


Spectrum 


18362009- 
2454424 



16 



Kunder et al. 



TABLE 3 

BRAVA Rotation and Dispersion Results 



gal / gal b Number Skewncss fr<,-^,^,„, Kurtosis (Jkuri p-indcx p-index 

(lU-g) ((1ck1 \'plori- " ' (Sliapiro- (Audprson- 

licH Wilkj Darling) 



-6.5 


-3.1 


63 


-0.60 


0, 


.31 


0.07 


0.62 


0, 


.043 


0.126 


5.5 


-3.1 


86 


0.05 


0, 


.26 


-0.31 


0.53 


0, 


.842 


0.991 


0.0 


-3.5 


111 


0.19 


0, 


.23 


-0.41 


0.46 


0, 


.035 


0.004 


4.0 


-3.5 


106 


0.22 


0, 


.24 


-0.64 


0.48 


0, 


.182 


0.266 


1.0 


-4.0 


61 


0.14 


0, 


.31 


0.86 


0.63 


0, 


.238 


0.267 


7.0 


-4.0 


93 


-0.08 


0, 


.25 


0.37 


0.51 


0, 


.386 


0.502 


6.0 


-4.0 


106 


-0.11 


0, 


.24 


0.09 


0.48 


0, 


.967 


0.943 


12.0 


-4.0 


104 


0.07 


0, 


.24 


-0.71 


0.48 


0, 


.203 


0.412 


14.0 


-4.0 


110 


0.28 


0, 


.23 


-0.23 


0.47 


0, 


.500 


0.306 


18.0 


-4.0 


103 


0.27 


0. 


.24 


-0.38 


0.48 


0, 


.434 


0.498 


22.0 


-4.0 


101 


-0.38 


0, 


.24 


0.72 


0, 


.49 


0, 


.258 


0.567 


-1.0 


-4.3 


109 


-0.17 


0, 


.23 


-0.09 


0, 


.47 


0, 


.513 


0.160 


4.0 


-4.5 


108 


-0.05 


0, 


.24 


0.02 


0, 


.47 


0, 


.750 


0.384 


6.0 


-4.8 


106 


0.26 


0, 


.24 


-0.37 


0, 


.48 


0, 


.357 


0.747 


-6.5 


-5.1 


79 


0.22 


0, 


.28 


-0.06 


0, 


.55 


0, 


.789 


0.728 


-10.0 


-6.0 


103 


0.15 


0, 


.24 


-0.36 


0, 


.48 


0, 


.440 


0.163 


-9.0 


-6.0 


102 


-0.02 


0, 


,24 


-0.50 


0, 


.49 


0, 


,899 


0.883 


-8.0 


-6.0 


105 


-0.24 


0, 


,24 


-0.26 





.48 


0, 


,449 


0.762 


-7.0 


-6.0 


103 


-0.04 


0, 


24 


1.96 





.48 


0, 


,087 


0.495 


-6.0 


-6.0 


96 


0.43 


0, 


25 


1.43 





.50 


0, 


116 


0.542 


-5.0 


-6.0 


108 


0.25 


0, 


24 


0.57 





,47 


0, 


,075 


0.329 


-4.0 


-6.0 


112 


0.38 


0, 


,23 


0.11 





,46 


0, 


,084 


0.426 


-3.0 


-6.0 


109 


0.01 


0, 


,23 


0.55 





.47 


0, 


,803 


0.506 


-2.0 


-6.0 


109 


0.06 


0, 


,23 


1.62 





.47 


0, 


,130 


0.364 


-1.0 


-6.0 


110 


-0.18 


0, 


,23 


0.31 





,47 


0, 


,898 


0.835 


0.0 


-6.0 


109 


0.09 


0, 


,23 


-0.62 





.47 


0, 


,360 


0.298 


1.0 


-6.0 


112 


0.22 


0, 


,23 


0.03 





.46 


0, 


,793 


0.614 


2.0 


-6.0 


110 


0.26 


0, 


,23 


-0.36 





.47 


0, 


,517 


0.754 


3.0 


-6.0 


114 


0.24 


0, 


,23 


-0.55 





.46 


0, 


,092 


0.109 


4.0 


-6.0 


113 


-0.33 


0, 


,23 


0.60 





.46 


0, 


,083 


0.556 


5.0 


-6.0 


109 


0.18 


0, 


,23 


-0.58 





.47 


0, 


,436 


0.523 


6.0 


-6.0 


102 


-0.35 


0, 


,24 


-0.09 





.49 


0, 


,104 


0.121 


7.0 


-6.0 


111 


-0.03 


0, 


,23 


-0.37 





.46 


0, 


,661 


0.507 


8.0 


-6.0 


107 


0.31 


0, 


,24 


0.14 





.47 


0, 


,572 


0.459 


9.0 


-6.0 


111 


-0.48 


0, 


,23 


0.03 





.46 


0, 


,680 


0.186 


10.0 


-6.0 


192 


-0.19 


0, 


,18 


0.79 





.35 


0, 


,153 


0.342 


-6.0 


-6.9 


97 


-0.12 


0, 


,25 


-0.80 





.50 


0, 


,121 


0.292 


6.0 


-6.9 


98 


0.32 


0, 


,25 


-0.22 





.49 


0, 


,133 


0.110 


-10.0 


-8.0 


97 


0.27 


0, 


,25 


0.07 





,50 


0, 


,622 


0.443 


-9.0 


-8.0 


104 


0.00 


0, 


,24 


-0.01 





.48 


0, 


,894 


0.864 


-8.0 


-8.0 


99 


0.24 


0, 


,25 


0.53 





.49 


0, 


,655 


0.545 


-7.0 


-8.0 


102 


0.41 


0, 


,24 


1.01 





.49 


0, 


,059 


0.117 


-6.0 


-8.0 


94 


1.38 


0, 


,25 


4.91 





.51 


0, 


,0001 


0.0001 


-5.0 


-8.0 


108 


-0.06 


0, 


,24 


-0.10 





.47 


0, 


,627 


0.649 


-4.0 


-8.0 


99 


0.78 


0, 


,25 


2.68 





.49 


0, 


,002 


0.227 


-3.0 


-8.0 


114 


0.01 


0, 


,23 


-0.32 





.46 


0, 


,466 


0.509 


-2.0 


-8.0 


102 


-0.42 


0, 


,24 


0.16 





.49 


0, 


,255 


0.320 


-1.0 


-8.0 


96 


0.71 


0, 


,25 


2.42 





.50 


0, 


,006 


0.079 


0.0 


-8.0 


102 


-0.01 


0, 


,24 


0.03 





.49 


0, 


,922 


0.773 


1.0 


-8.0 


103 


-0.05 


0, 


,24 


0.40 





.48 


0, 


,762 


0.686 


2.0 


-8.0 


103 


0.21 


0, 


,24 


-0.09 





.48 


0, 


,467 


0.268 


3.0 


-8.0 


108 


-0.33 


0, 


24 


0.89 





,47 


0, 


,422 


0.299 


4.0 


-8.0 


104 


0.11 


0, 


24 


-0.13 





,48 


0, 


,802 


0.934 


5.0 


-8.0 


102 


0.24 


0, 


,24 


-0.05 





.49 


0, 


,519 


0.481 


6.0 


-8.0 


98 


0.43 


0, 


,25 


0.29 





.49 


0, 


,173 


0.110 


7.0 


-8.0 


102 


-0.23 


0, 


,24 


0.60 





.49 


0, 


169 


0.182 


8.0 


-8.0 


99 


-0.89 


0, 


,25 


4.16 





.49 


0, 


,0001 


0.090 


9.0 


-8.0 


108 


-0.61 


0, 


,24 


1.89 


0, 


.47 


0, 


,007 


0.288 


10.0 


-8.0 


96 


0.06 


0, 


.25 


-0.61 


0, 


.50 


0, 


.460 


0.585 



